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Abstract 

Based on a quasiparticle model for /3 stable and electrically neutral deconfined matter we address 
the mass-radius relation of pure quark stars. The model is adjusted to recent hot lattice QCD 
results for 2 + 1 flavors with almost physical quark masses. We find rather small radii and masses 
of equilibrium configurations composed of cold deconfined matter, well distinguished from neutron 
or hybrid stars. 
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I. INTRODUCTION 



After growing evidence for the quark-gluon substructure of hadrons the question has 



been asked jl whether massive neutron stars may have a core composed of quarks 



These so-called hybrid stars may be part of the neutron star branch or constitute a separate 
stable branch of high-density objects - the so-called third family S| or twin stars flo| . 
Also pure quark stars populating another separate branch of stable, spherically symmetric 
cold objects have been discussed , Il2j. All these possibilities depend sensitively on 



the equation of state at high density and the details of the deconfinement transition at low 
temperature. While at high temperature and zero net baryon density a proper numerical 
evaluation of the equation state based on first-principles - QCD - is accomplished, the 
knowledge of the equation of state at high baryon density and low temperature is fairly 
poor. In the asymptotic region, safe statements on the matter states can be made |l3, Q, 
but the extrapolation to the interesting region of energy densities around 10 15 g/cm 3 is 
hampered by serious uncertainties as one expects significant non-perturbative effects. 

A possibility to approach the theoretical analysis of quark stars is to employ certain mod- 
els adjusted to high-temperature lattice QCD results at zero or small net baryon density. 
Of course, the applicability of such models at low temperatures and high densities is not 
guaranteed. Quarks stars or neutron stars with quark cores are expected to have similar 
mass-radius relations as ordinary neutron stars. This makes difficult an experimental veri- 
fication via these observables. The modified cooling behavior of quark matter is considered 



15|. 



as a possible tool to find appropriate observational hints 

Here we are interested in the mass-radius relation of pure quark stars which are cold and 
spherically symmetric. We rely on a quasiparticle model (cf. [I6MI8I] for such models) which 
we adjust to recent realistic lattice QCD results. Our quasiparticle model |l6,Q, Q| allows 
for a suitable parametrization of lattice QCD data at zero and non-zero chemical potential. 
Its structure can be derived from a two-loop $ functional 21 d3|- To accommodate further 
non-perturbative effects the running coupling g s is replaced by an effective coupling G. In 
the simplest version the imaginary parts of the self-energies are neglected and the dispersion 
relation is approximated by utilizing the asymptotic self-energy. The model has been shown 
to describe successfully various lattice QCD data at zero chemical potential, at non-zero 
(including also purely imaginary) chemical potential of bulk thermodynamical quantities up 
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to off-diagonal susceptibilities 



24- 
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28fl are 



New high-temperature lattice QCD data for almost physical quark masses 
now at our disposal for zero chemical potential. We adjust our model at this data and 
extrapolate the equation of state to zero temperature. The emerging equation of state is 
;hen used to consider cold pure quark stars. Analog studies have been performed in, e.g., 



29H33J, however without such intimate contact to advanced lattice QCD results. 

Our paper is organized as follows. In section [IT] we formulate our model for zero tem- 
perature. The comparison with hot lattice QCD results is performed in section IHH The 
parameters are used in section [TV] to gain the cold equation of state. The emerging mass- 
radius relations of cold equilibrium configurations are discussed in section [V] The summary 
can be found in section [VTl The Appendix lists expressions used for transferring the hot 
lattice QCD data to finite baryon densities. 

II. QUASIPARTICLE MODEL AT T = 

For the employed quasiparticle model the pressure p = J2 i=u ds Pi and quark densities rii 
at temperature T = are given by 



PM = ^ / dk : - Bi{fMi), (1) 

to* Jo ^/k 2 + m 2 

B^) = Bi(no) + / dp. — ^ / dk ■= (2) 



*x 2 Jn 9fl J ^k 2 + m 2 ' 

n^) = ^-M~ m "f /2 (3) 

with the index i denoting the quarks u, d and s with degeneracies d u ^ s = 2iV c = 6. The 
choice of /i and the corresponding integration constant B(/jl ) is described below. The 
energy density follows from e = J2i= u ds(Pi + Hi n i)- The asymptotic quark masses, which 
enter the employed dispersion relations ojf = k 2 + m 2 as approximation of the self-energies, 
are 

ml = m 2 i)0 + 2m ifi Mi + (4) 



M 2 = ^[T 2 + 4)G 2 (5) 



n 2 
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with Cf = (N 2 — l)/(2N c ), where the rest masses m it0 may be included accordingly. We 
employ m Uj o = map = rn S fi/10 with m S) o = 105 MeV as in [27|, |28fl. For later use also the 
temperature dependence is displayed here and we note already the gluon g asymptotic mass 
at finite temperature 




i=u,d,s / 



+ ^ Y, ^]G 2 , (6) 



where = N c + | and the degeneracy factor d g = N 2 — 1. 
The five relations for charge neutrality 
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^n u (fi u ) - -n d (/j, d ) - -n s (fi s ) - n e (fi e ) - ra M (jU M ) = 0, (7) 

for P equilibrium 

Hd = l^u + He (8) 

(e.g., from n «-> p + + e~ + u e ), for equilibrium due to strangeness changing weak decays 



(e.g., from A <-> p + + n ), for \i decay 



/is = /id (9) 



= /i e (10) 



(e.g., from /i «-> e + v e + z/^), and for total baryon density 

n(/i) = ~(n u (// w ) + n d (/i d ) + n, (//«)) (11) 

map the various chemical potentials on one independent baryon chemical potential /i via 
/iu,d,s,e,^(/i) as required as consistency condition of the utilized quasiparticle model. We 
assume that the neutrinos v e [1 left the star matter and, therefore, do not participate in the 
chemical equilibrium reactions. The pressure and density expressions for the electron e and 
muon \i components are as Eqs. ([TJ) (without the functions Bj) and ((3l). 



The effective coupling G 2 follows from the flow equation 



0,0] 



dG\^T) dG 2 (fi,T) 

aT or +a » dli - 6 ' (12) 

where the coefficients a-r, and b (cf. Appendix lAl) depend again on the effective coupling 
G 2 as well as both temperature and chemical potential. It is integrated using the method 
of characteristics. Along each characteristic line, the input information G 2 (T), extracted 



from lattice QCD data in the next section, is transported from the temperature axis to the 
chemical potential axis thus providing G 2 (fi). Along one arbitrary characteristic, emerging 
at T = T and meeting the /i axis at fio, the meanfield contribution B is integrated (as 
outlined in Appendix lAl too) . yielding the necessary integration constant B(ft ). With the 
effective coupling G 2 (fi) and the B(fio), all thermodynamic quantities along T = are then 
determined. 



III. EQUATION OF STATE FROM LATTICE QCD DATA AT /x = 



In [28j] ([2 711) the interaction measure A(T)/T 4 = (e — 3p)/T 4 has been presented for 
almost physical quark masses for the two light quarks and a strange quark in the temperature 
range T = 140 - 475 MeV (140 - 825 MeV) at fn = 0. We rely here on the p4 and asqtad 
data for N T = 8 1281 and 6 [271 and assume that further cut-off effects are negligible, i.e. we 
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compare our continuum model with the finite-size results [2T|, |28j| . It seems most appropriate 
to adjust our parameters directly at the interaction measure, being the primary information 
from lattice QCD, which reads in the quasiparticle model [19( for fi — 

v (is i rw (* 2 + m .vM (13) 



2tt 2 J To dT' J eV /fc2+ mi (T') 2 /r' + 5 . 

where S g = —1 and S Ut d, s = 1 and mj from Eqs. (jlj) and (jBJ) with explicit and implicit T 
dependencies. The latter one is in G 2 for which we choose 

with ^ = T ~ Ta as well as (3$ = 11 — liVf for iV/ = 2 + 1 flavors as a convenient parametrization 
of the effective coupling which resembles a regularized 1-loop running coupling. 

Using general thermodynamic relations one can calculate the pressure via p(T)/T 4 = 
p(T )/Tq + J^ o dT' A(T') T'~ 5 , where p(T ) is an integration constant. The chosen value of 
T should be at the lower limit of our model for deconfinement, i.e. T ~ 190 MeV according 



to [28J]. The scaled entropy density is accordingly s/T 3 = 4p/T 4 + A/T 4 ; unfortunately, it 
depends also on the pressure normalization via 4p(T )/T 4 . 

Our quasiparticle model is based primarily on the entropy density s, i.e. pressure and 
interaction measure are analytic integrals of the entropy including an integration constant 



action N T T s [MeV] A [MeV] 


— Fir, 


v 2 /rlnf 


p4 6 


167 


20 


-(78 MeV) 4 


0.821 


p4 8 


146 


31 


(163 MeV) 4 


0.679 


asqtad 6 


131 


45 


(90 MeV) 4 


0.981 *) 


asqtad 8 


107 


61 


(166 MeV) 4 


0.654 



Table I. Parameters of fits to the lattice data 
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28|] for T > 190 MeV. B = £\ = 



-g,u,d,s Bi,0- 



*) Without the data point at 213 MeV which would drive the fit to fail in the high-temperature 
region. 

B = ^2Bi tQ . Thus fitting the interaction measure means not only determining the pa- 
rameters T s and A of the effective coupling but also the pressure integration constant B . 
Thus, our pressure as well as entropy and energy density follow directly from e — 3p without 
another integration constant. (This is due to additional knowledge of explicit expressions 
for all thermodynamic quantities, as opposed to general thermodynamic relations, where an 
additional constant p{T ) is required to arrive from the interaction measure at the pressure. 
Within the quasiparticle model, p(T ) is known from the parameters T s , A and B Q via the 
expressions of s and A in p(T Q ) = (Ts(T ; T s , A) — A(T ; T, A, £?n))/4.) 
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2811 to the interaction measure 



The x 2 minimization of the difference of the data in 
Eq. (fl~3l directly yields the values of T s , A and B as listed in Tab. fl] With this parametriza- 
tion we get the interaction measure as exhibited in the left panel of Fig. [TJ The maximum 
of A arises from a turning point of the scaled pressure as a function of log T. Within our 
quasiparticle model, the location of the maximum is governed by the values of T s and A, 
where the latter one also affects the peak width. The peak height of A/T 4 is essentially de- 
termined by Bq. The fits are in a narrow corridor for T > 300 MeV yielding some confidence 
in the equation of state there. 

of Fig. Q] we compare the pressure of our model with the pressure 



In the right pane 
estimate deduced in 

heights in A/T 4 , the resulting pressures in our model are in a reasonably narrow corridor: 



281 ] from the interaction measure. Despite of the variation of the peak 



Our fit to the asqtad N T = 8 data is in the middle of the pressure range determined in [28| by 
different interpolations on the data for A/T 4 and assumptions on po. The p4 N T = 6 peak 
in A/T 4 is higher, and, consequently, our pressure is also somewhat higher, governed by the 
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Figure 1. (color online) Left panel: Fits of the quasiparticle model interaction measure Eq. ([13]) to 
the lattice QCD data (symbols) for lattice actions p4 (blue) and asqtad (red) and lattice spacings 



N T from 
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281 ] . Solid (dashed) curves are for N T = 8 (6). Right panel: Scaled pressure p/T 



of the quasiparticle model adjusted in the left panel compared to the pressure estimate from 
(grey area). 



28] 



positive B . A fit to the upper and lower limits of the pressure band from 



28J] would result 



in (T S ,X,-B ) = (109 MeV, 53 MeV,(185 MeV) 4 ) and (36 MeV, 107 MeV, (197 MeV) 4 ), 
respectively. 



IV. THE QUARK EQUATION OF STATE AT ZERO TEMPERATURE 

Utilizing the values of Tab. [J the flow equation lfi~2l) is solved using the mentioned method 
of characteristics to obtain G 2 (fi) and B(fi). In doing so, the side conditions (l7T fTTl) are 
invoked so that along each characteristic curve the requirements of ft stability and electric 
charge neutrality are fulfilled. The characteristics for the p4 action and N T = 8 are exhibited 
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in Fig. El As already noted in [12j, |l_9|] , the characteristics emerging from the very vicinity 
of T have the tendency to cross each other at low temperatures. (An extended version of 
the model in cures this insanity.) The pressure becomes negative at [i < 550 MeV. 

Clearly, we consider only the region of positive pressure where the characteristics behave 
regularly. 

It happens that the effective coupling G 2 can also be parametrized at vanishing tempera- 



7 





200 400 600 800 

H [MeV] 



Figure 2. Characteristics of the flow equation (1121) with side conditions (|7lfTT1) imposed for the p4 
action and iV r = 8. 

ture using Eq. (fT5l ) but with £ = The parameters are listed in Tab. HT1 The interaction 

measure A(//)/// 4 displays a peak, as A(T)/T 4 does. 

The pressure contributions according to Eq. {TJ with the such obtained effective coupling 
G 2 (/i) are exhibited in the left panel of Fig. [3] as a function of the chemical potential /i = ji u . 
The differences of up, down and strange quark contributions are determined by differences in 
the respective chemical potentials as shown in the right panel of Fig. [31 Due to equilibrium 
with respect to strangeness changing weak decays, n<i = Us holds which deviates slightly 



from fi u . In line with [IJ] the lepton contributions are tiny, as evidenced in Fig. [3], too. The 
pressure difference of down and strange quarks is due to the considerably larger rest mass 
of the latter ones. For the sake of completeness we also show the individual contributions 
to energy density (left panel in Fig. BJ and the individual particle densities (right panel in 
Fig. H|). Similar to the pressure, the lepton contributions are not visible on the used scales. 

The resulting equation of state at T = in the form e(p), needed for the integration 
of the TOV equations below, is exhibited in Fig. [5] (left panel for a comparison of the four 
equations of state) together with two fits by e = v~ 2 p + e for the equations of state adjusted 
to N T = 6 and 8 for the p4 action (right panel). Parameters for all actions and temporal 
lattice extends considered here are listed in Tab. IHIl Both, the vacuum energy density 
e o = e(p = 0) and the velocity of sound parameter v 2 = dp/de are in narrow intervals for 
the fours sets of lattice QCD input data: While the vacuum energy density varies within 
(366 MeV) 4 - (381 MeV) 4 , v~ 2 is within 3.8 - 4.5. As the interaction measure A(/i)//i 4 
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action N T fj, s [MeV] A M [MeV] 



p4 6 211 159 

p4 8 134 215 

asqtad 6 68 285 

asqtad 8 -72 380 

Table II. Parameters of Eq. (fT5l) with £ = M ^ Ms following from the solution of the flow equation 



Eq. (H2|). The fits apply in the range y = 0.6... 1.2 GeV. 
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Figure 3. Left panel: The scaled pressure contributions Pi/y 4 for the p4 action and N T = 8 as 
functions of /i. The leptonic contributions are on the lower /j, axis in the given Pi//J, scale. Right 
panel: The individual chemical potentials fii as functions of // = fj, u . 

becomes small at large values of /i, our resulting equations of state e{p) have the tendency 
to merge. (Some differences are caused by the different fit values of -Bo-) At small pressure 
the deviation of our four equations of state are about 10%. 

The same is true if considering again the upper and lower limits of the pressure band 



from 



28| instead of the interaction measure. We find values (v s 2 , e 



1/4^ 



(3.92, 370 MeV) 



and (4.24, 383 MeV) in the parameter area of the above fits where larger values of the 
pressure at vanishing chemical potential lead to smaller values of the vacuum energy density 
at vanishing temperature and the inverse squared velocity of sound. Also, thermal effects 
are found to be small, i.e. up to T = 50 MeV the equation of state e(p) does not change 
significantly. 

Our results can be compared with perturbative calculations at vanishing temperature in 
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Figure 4. The scaled energy density contributions ei/fj, 4 (left panel) and the scaled net quark density 
rii/fx 3 (right panel) for the p4 action and N T = 8 as functions of //. For the used scale, the leptonic 
contributions are not visible. 

action N T v~ 2 e^ 4 [MeV] 
p4 6 3.81 381 
p4 8 4.01 366 
asqtad 6 4.23 379 
asqtad 8 4.47 367 

Table III. Parameters of linear fits e = v~ 2 p + eo to our equation of state with G 2 (/j.) determined 
by the flow equation (fl~2l) . The leptonic contributions are included. 
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321 ] . In [32|] the pressure of cold quark matter is calculated in hard-dense-loop pertur- 



bation theory. The resulting pressure for 3 flavors with equal chemical potential and the 



choice [32|] of the renormalization scale p = p is shown in the left panel of Fig. [6l Also 
shown are the results from the weak-coupling expansion to second order fsi]. Here the value 



of p is varied from p to 2p. For reference also a comparison with NJL model results 
is depicted. The grey dotted curve is p = 3 GeV/fm 3 ; the region < p < 3 GeV/fm 3 is 
relevant for quark stars, as turns out by integrating the TOV equations (see section IVT). 

The energy density e = ji 2 d/djj,{p/ ' p) depends on the incline of the pressure scaled with 
the chemical potential rather than the absolute values of the pressure. This explains the 
fact that, while the scaled quasiparticle pressure p/p 4 from lattice QCD (left panel in Fig. [6]) 
shows some spread as a function of p, the resulting equation of state e(p) (right panel in 
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Figure 5. (color online) Left panel: The energy density as a function of the pressure at T = 



(solid curves) for all four considered lattice results from 
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281 ] . Blue (red) lines represent results 



deduced from p4 (asqtad) lattice data while solid (dashed) curves denote N T = 8 (6). Right panel: 
Two of the linear fits e = v~ 2 p + e$ (dotted curves) from Tab. II III compared to the corresponding 
equations of state (color code as in left panel). 



Fig. [6]) is given as a tight band. Inspection of p//i as a function of \i (not displayed) explains 
the broad range of values for e(p — > 0) in [311 ]: the slope of as a function of // changes 
drastically with the chosen scale Jx for small pressures. For // = 1.5// the equation of state in 



31] in the form e(p) coincides with the results of [32], which in turn falls in the same range 
as our set of equations of state. In fact, v~ 2 = 3 and ejj/ 4 = 365 MeV yield a good description 
of the equation of state from JsjJ for // = 1.2// and 32]. (Expanding the quasiparticle partial 
pressure pi (pp) including the meanfield contribution Bi ([2]) at T = in powers of the coupling 
constant G yields the leading terms Pj(//j) = (1 — 2a s /n + . . .) (47r 2 ) + £?j(/i ), where the 



coefficient of the 0(a s ) term, a s = 4nG 2 , equals the strictly perturbative results in [3lL l35| 



the coefficient of the next-order term deviates from t he p erturbation expansion, similar to 



the quasiparticle model 
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hard-dense-loop approach in 



22l.) 



161] ( see also discussion in |22|]) at non-zero temperature and the 



Remarkable is that all the discussed equations of state have a certain value of the chemical 
potential at vanishing pressure. This enables, in principle, to construct pure quark stars with 
vanishing pressure at the surface. 

Fig. [6] clearly evidences that the previous foundation for discussing quark stars seemed 
not to be on safe grounds as the proposed model equations of state were too different 
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Figure 6. (color online) The scaled pressure p/fi 4 as function of n at T = 0( 



equation of state e(p) (right panel) for the four considered lattice results from 
32] (dash-dotted curve), 33]] (dash-double-dotted curve) and 



to results from 



27 



eft panel) and the 



281 ] in comparison 



311 ] (grey bands limited 



by fi/n = 1, 1.5 and 2 from dark to light). Blue (red) curves represent quasiparticle results adjusted 
to p4 (asqtad) lattice data while solid (dashed) curves denote N T = 8 (6). Leptonic contributions 
are included. 

unless further constraints (as the compatibility, e.g., with a hadronic model equation of state 



required in 



3l| ) are imposed. Given the intimate contact of our approach to first-principle 



evaluations of QCD, we hope to have a more reliable foundation. Of course, this hope is 
related to the assumption that the extrapolation to non-zero chemical potential is sufficiently 
smooth. The successful comparison of our model with Taylor expansion coefficients for the 



\x dependence 



24J as well as the application of our model at imaginary chemical potential 



251 ] (not only small values thereof!) give us some confidence in our approach. 

Let us finally comment on the importance of the side conditions. If one assumes one 
common chemical potential for all quarks [i and includes leptons (// e = /i M ) via a electric 
neutrality condition /i e = /x e (/-0 , the results of the equation of state differ from the isospin 
asymmetric model with the side conditions (l7lfTT!) properly invoked on a 10% level. 
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V. INTEGRATION OF THE TOV EQUATIONS 



To estimate the properties of quark stars as spherical equilibrium configurations of pure, 
strongly interacting quark matter we employ the TOV equations 

dp _ ([1 + v~ 2 ]p + e ) (m + 4vrr 3 p) 

dr r 2 ( 1 - ^fG N ) 

^ = ^r 2 (v; 2 p + e ), (17) 

where the special parametrization e = v~ 2 p + e Q of the equation of state is supposed to hold. 
Gat is the Newtonian gravitational constant, and we employ units with he = 1. 

We emphasize the strong dependence on the actual value of eo which determines the 
pressure gradient in the dimensionless combination Gnc]^ 2 (which is of the order of 10~ 39 
for the case at hand), which can be seen in writing the TOV equations as 

dp ([1 + v~ 2 ]p + 1) (m + 47rf 3 p) 

dr f 2 (l-fO ' 1 ' 

— = 47rf 2 (v; 2 p + l), (19) 

from the scaled quantities p = peo, r = f(GAreo) -1 ' 2 , m = m (GWeo) -1 / 2 ^ 1 . The scaled 
TOV equations depend only on v~ 2 . The solutions for the relevant values of v~ 2 = 2. ..4 are 



exhibited in Fig. [71 With the given scaling, m and r shrink with increasing value of e 



1/2 

5 



while the dependence on v~ is moderate within the interval covering the values of Tab. | 
Thus the vacuum energy density eo is indeed the decisive quantity determining the sizes and 
the masses of pure quark stars. To be specific, for v~ 2 = 3 ± 1, the scaled maximum mass 
is 0.004 ± 0.001. 

To test the dependence of deviations from the approximation e = v ~ 2 p + e we integrate 
the TOV equations with our equations of state adjusted to the lattice QCD results. The 
results are exhibited in Fig. [3 The maximum masses are about O.5M with radii of about 3 
km. If such objects would exist, their bulk characteristics were quite different from canonical 
neutron stars with masses concentrated at 1.4 M and radii of 15 km and larger. Therefore, 
;he pure quark stars from our analysis cannot serve as candidates of twin stars discussed in 

We stress again the important role of the value of eo = e(p = 0). With the above derived 
scaling, equations of state with significantly smaller values of e than deduced in our analysis 
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of the lattice QCD results combined with the employed quasiparticle model, would allow for 
significantly larger masses and radii. 

The present considerations will be modified when combining our equation of state of 
deconfined matter with a hadronic low-density equation of state at p > 0. Then hybrid stars 
could be constructed with properties depending to a large extent on the transition region 
from confined to deconfined matter. 
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Figure 7. Scaled mass M = m(p = 0) shown as a function of scaled radius R = r{p = 0) for several 



values of v s as solution of the scaled TOV equations (|T8l) and (1191) . 



VI. SUMMARY 

In summary we employ a quasiparticle model, adjusted to recent realistic lattice QCD 
data with almost physical quarks masses, to consider pure quark stars. The needed equation 
of state can be approximated very well by the concise form e = v ~ 2 p + en with values of 



28j. Lattice data 



u- 2 = 3.8 - 4.5 and ej /4 = 366 - 381 MeV from the lattice QCD data 
from both the p4 and the asqtad version can be described equally well and lead to similar 
spherically symmetric stars. The maximum masses are about 0.5 M with radii of 3 km. 

The pure quarks star masses and radii scale with e$ which is the decisive quantity as 
eo is the vacuum energy density at vanishing pressure. It follows within our model directly 
from the lattice QCD data at finite temperatures. 

Rapidly rotating quark stars exhibit a disc like shape with sharp edge. Their maximum 
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results 
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or online) Mass-radius relations of quark stars following directly from lattice QCD 
281 ] within our quasiparticle approach. Blue (red) lines represent results deduced from 
p4 (asqtad) lattice data while solid (dashed) curves denote N T = 8 (6). 

at the shedding limit. 



20 



231 ] . where 



masses are enlarged by 65% and the radii by a similar amount 

The present approach can be extended to the full HTL quasiparticle model 
effects of Landau damping and collective modes are included 
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Appendix A: Coefficients of the flow equation 

The coefficients in Eq. ( fi~2"l) are 



(It 



- E 

i=u,d,s 

E 

i=g,u,d,s 

E 

i=g,u,d,. 



drii dm 2 








dm 2 dG 2 

s 1 








dsi dm 2 








dm 2 dG 2 


j 






drii dm 2 




dsi 


dm 2 


dm 2 dT 


G 2 ,m 


dm 2 


<9/i 



G 2 ,T 



with 
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drij 
dm 2 
dsj 
dm 2 



di 



4vr 2 T J 
di 



An 2 T 



OO 7,2 

_ _ IX „ 

dk—F. 



dk k 2 



where F± = [f 2 exp ^'"^ ± /£ exp and /± are the statistical distribution functions 

for fermions (+) and anti-fermions (-) respectively. The derivatives of the effective gluon 
masses (l6l) are 



dm 2 g 



dT 

dm 2 g 



G 



dfx 



G 2 ,T 



127T 2 



2 (Ni+N h ) n + (Ni+2N h ) HejG 2 + (Ni + 2N h ) (// + // 



<9/i e 



9/i 



G 2 ,T, 



9m 2 



«9G 2 



brp2 



T + 
6 127T 2 



(iVi + 2JV fc ) (fi + fi e 



dfi e 



dG 2 



G 



where N\ and N h are the numbers of included light (2) and heavier (1) quark flavors. For 
the effective quark masses (i = u, d, s) one has 



dm 2 




dm 2 


+ 


dm 2 




dfi e 




dT 




dT 




dfi e 


G 2 ,T,/i 


dT 


G 2 ,)M 


dm 2 




dm 2 


+ 


dm 2 




dfi e 




dfi 


G 2 ,T 


dfi 


G 2 ,7> e 


dfi e 


G 2 ,T,fi 


dfi 


G 2 ,T 


dm 2 




dm 2 


+ 


dm 2 




dfie 




dG 2 


T,fM 


dG 2 




dfi e 


G 2 ,T,im 


dG 2 


2> 



with dm 2 /dT\ G2 
and dmj/dG 2 



VTG 2 , dm 2 /dfi\ G2T = V/aG 2 /* 2 , dm 2 /dfi 



e\G 2 ,T,n 



5 iu V fiiG 2 / H 2 



| T ^ e = (T 2 + fit/n 2 ) /2, where V = (m q>0 /M q + 2) C f /4. The derivatives of 
the electron chemical potential therein are 

dfi e 

— — vv 

G 2 ,A* j=u,d,s,e,fj, 



dT 

dfie 



dn,; 



dT 



G 2 



dfi 
dfi, 



G 2 ,T 



drij 



j=u,d,s,e,/i 



d/i 



dG 2 



W v q , 9n ' 



dG 2 



G 2 ,T,ii e 
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with W = Ej=M,»,e,^* dn j/ d ^\ G 2^ T and 



dT 

drij 
d\x 

drii 



G 2 , M, Me 



G 2 ,T. 



/'<• 



<9/i e 

<9m 



<9G 2 



G 2 ,fi,T 



T,fi,fj, e 



dT 



drij 
dm dm 2 



+ 



drij dmf 



... am 2 <9T 
9n,: (9m 2 



G 2 ,M, /ie 



G2jm 2 <9m 2 

9n,- (9m 2 



G 2 ,T. 



Me , 



<9/X 



+ 



G 2 ,m 2 



dmf d\ii 



G 2 ,/i. 



d[M e 



dmf dG 2 



as well as 



duj 
dT 

dn 
(9/i 



di 



G 2 ,fi,iie,m 2 



G 2 ,m 2 



, 2tt 2 T 



poo 

/ d£; A; 2 - HiF+} 

Jo 



di 



2vr 2 T 



dk k 2 F+. 



qi are the electric charges of the quark species. Note that the side conditions Eqs. (l7lfITl) are 
included. These strongly modify the coefficients given in Jj]]. 

Along the characteristics, where /x, T and G 2 are given as functions of the affine curve 
parameter x, the bag pressure B has to be integrated according to 



B = B(n = 0)-Y, r 
,• Jo 



dx 



dpi 
dm 2 



dm 2 



or 



dT 



dm 2 



g 2 ,m 



+ b 



dm 2 



G 2 ,T 



dG 2 



with 



dpi 
dm 2 



di 



47T 2 



poo 7,2 

/ dk-[f + + f-]. 
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